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I. INTRODUCTION 

During the last decade two-proton emitters and 
proton-rich nuclei, in the vicinity of the proton drip- 
line, are subjects of intensive experimental and theoret- 
ical studies. In particular the study of a candidate to 
possess a two-proton halo in the ground state, namely 
the Borromean 17 Ne nucleus is of special experimental 
and theoretical interest (see, for example, discussions in 
[E IE 0)- The inherent feature of the halo structure is 
a relatively small separation energy of a valence nucleon. 
It reveals itself in a large valence nucleon removal cross 
section and a narrow core longitudinal momentum (LM) 
distribution. 

In the case of 17 Ne the proton removal cross section, 
measured at the energy 66 A MeV on a Be target |2(, is 
relatively large compared to the cluster model ( 15 0+p-|- 
p) predictions |3j. At the same time, the measured core 
LM distribution is wider than the calculated one. Both 
these facts can be attributed to contribution of a proton 
removal from the 15 core in 17 Ne if this cross section is 
relatively large. 

Recently the 14 N longitudinal momentum (LM) distri- 
bution and break-up cross section (into the 14 N+p chan- 
nel) have been measured in fragmentation of 15 O on a 
Be target at the energy 56 A MeV 0. This opens a 
possibility for more precise calculations of the proton re- 
moval from 17 Ne and evaluations of the contribution of 
the proton removal from the 15 O core to this process. 

In this paper, we present a detailed analysis of the 15 O 
break-up in light targets. We perform the calculations in 
the eikonal approximation of the Glauber model |5j, |6J, |7j, 
8]. This approach is well developed and convenient for 
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calculations of break-up cross sections, interaction cross 
sections, and momentum distributions of fragments in 
break-up of a nucleus at intermediate and high energies 
(from 30 to 1000 A MeV). 

The formalism for the calculations is described in Sec- 
tion [HJ The main ingredients of the Glauber model are 
the wave function of the relative motion of the fragments 
and the profile functions defining the fragment-target in- 
teraction. They are fixed using experimental data on 
the nucleon-nucleus and nucleus-nucleus cross sections, 
proton separation energies, the level scheme of the core 
nucleus, etc. In particular, the profile functions are fit- 
ted using the nucleus-nucleus and nucleon-nucleus inter- 
action cross sections. 

The wave function is obtained in the core+proton 
( 14 N+p) model of 15 O, where the 14 N core fragment can 
be in the ground and excited states (see, for example, Ref. 
@)- The p-wave proton removal from 15 (J 71 " = l/2~) 
leads to few 14 N bound states. We consider four of them 
(for details see @): E x = 0.0 MeV (J* = 1+, T = 0), 
E x = 2.313 MeV (0+, 1), E x = 3.948 MeV (1+ 0), and 
E x = 7.029 MeV (2+, 0). Here, E x is the excitation 
energy and (J 7r , T) are the spin and isospin of the 14 N 
state. For each state, the depth of the ( 14 N+p) inter- 
action potential (see, below) is fitted to reproduce the 
proton separation energy. 

The cross sections of the proton removal from the 15 O 
ground state are determined by the spectroscopic factors 
[3, El °f the p-wave proton states. 

In Section IIIII we fit the profile functions in calcula- 
tions of the corresponding nucleus-nucleus and proton- 
nucleus interaction cross sections, and compare results 
to the available experimental data. 

In Section IIVI we present the calculated cross sec- 
tions and longitudinal momentum distributions of the 
14 N fragments produced in various states in the process of 
the one-proton removal from 15 on a 9 Be target. These 
results are compared to the experimental data on the 
14 N longitudinal momentum distribution and the break- 
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up cross section measured at the energy 56 A MeV 



II. CROSS SECTIONS AND MOMENTUM 
DISTRIBUTIONS 

In the core-nucleon model of the projectile nucleus, 
the initial state is described by the wave function (WF) 
^ JMj (f) of the core-nucleon relative motion with a total 
angular momentum J and its projection Mj. The WF 
depends on the relative coordinate r between nucleon and 
core. 

After interaction with a target, the WF of the projec- 
tile will be corrected by factors, connected with nucleon- 
target and core-target interactions. Thus, the WF in the 
projectile rest frame is modified as Q 



where R is the coordinate of the center of mass of the 
projectile, bi — \bi\ (i = n,c), and b n , b c are the trans- 
verse two dimensional impact parameters of the nucleon 
and the core with respect to the target nucleus, i.e. 
b n = R± + f ± A c /(A c + 1) and b c = R ± - f±/(A c + 1), 
where R± and r± are components, perpendicular to the 
beam direction taken as z axis, and A c is the mass num- 
ber of the core. The profile functions S n (b n ) and S c (b c ) 
are generated by nucleon and core interactions with the 
target nucleus. 



V(f,R) = S c (b c )S n (b n )^ JMj (r), 



(1) 



The fragmentation includes nucleon stripping and 
diffraction processes. The corresponding cross sections 
are given by the equations [f| 
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The proton removal cross section is found as the sum 

C-p = (T s tr + &diff ©• 

The wave function ^ jm, is 



^JMj = [^LM{r)®Xs 



Xs 



JMj 



(3) 



where Xs c m c is the internal wave function of the core 
including the spin function, Xs n m n is the spin function of 
the valence nucleon. 

We denote the part of the WF related to the relative 
motion as \&lm 



^lm(t) = R L (r)Y LM , 



(4) 



where Ylm l is the spherical function. 

The radial part of the core-proton WF, Rl(t), is ob- 
tained as a solution of the Schrodinger equation for the 
Woods-Saxon potential (the Coulomb 14 iV + p potential 
is also included). For each state of 14 N, the parameter 
Vq of the Woods-Saxon potential is fitted to reproduce 
the proton separation energy with the fixed parameters 
a = 0.65 fm and R = l^A^^S.OO fm. The depth 
parameters and the proton separation energies are given 
in Table IJ 

In the calculations of the cross sections and LM dis- 
tributions of the fragments we consider the p (pi/2 and 



P3/2) proton removal. The p-wave proton removal from 
15 leads to the residual 14 N core in the bound states 
E x = 0.0 MeV (J* = 1+, T = 0), E x = 2.313 MeV (0+, 
1) E x = 3.948 MeV (1+, 0), and E x = 7.029 MeV (2+, 
0). 

Note, that the spectroscopic factors are not measured 
yet for 15 0. As was shown in the DWBA analysis the 
contribution of the protons with I = 1 dominates in the 
proton transfer reaction leading to the ground state of 
15 O. The spectroscopic factors of the states can be taken 
as those predicted by Cohen and Kurath [l(J . These val- 
ues are close to the measured values in the neutron pickup 
reactions with the mirror 15 N nucleus, the 15 N(p, d) 14 N 
reaction with 40 MeV protons |j , a nd the 15 N(d,£) 14 N 
reaction with 90 MeV deuterons [13] (see the discussion 
in Q and references therein). We use the spectroscopic 
factors from 01 an d E3- These factors C 2 S are also 
listed in Table [Q 

Note, that the contribution of the 14 N excited bound 
states to the diffraction cross section (j2J) is relatively 
small and is neglected here. 

The LM distributions of the core fragments are ob- 
tained by the Fourier transformation of the core-proton 
WF, i?^(r), corrected for the core-target and nucleon- 
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TABLE I: The depth parameter Vq of the Woods-Saxon po- 
tential, obtained with the diffuseness parameter do = 0.65 fm 
and radius Ro ~ 3.00 fm for the p-wave proton separation en- 
ergy E s . E x is the corresponding 14 N core excitation energy. 
C 2 S are the spectroscopic factors, a [TH ] and 6 (T!3 ]. 











Woods-Saxon potential 


E x 


"N 


C 2 3 a 


C 2 S b 


V 


E s 


MeV 


(J*,T) 






(MeV) 


(MeV) 





(l + ,0) 


1.459 


1.343 


-48.09 


7.297 


2.313 


(o + ,i) 


0.418 


0.472 


-52.07 


9.610 


3.948 


(l+,0) 


0.696 


0.656 


-54.78 


11.245 


7.029 


(2+,0) 


1.250 


1.250 


-59.73 


14.326 



target interactions 
da s tr 1 



dk z 2L + 1 



b n db n {\ - \S n {b n )\ 2 ) 



r ± dr ± dcj)\S c (\b n - fj_|)| 



(5) 



E 



J e tk ^R L (y/rl+z*)Y LML dz 



The core longitudinal momentum distribution in the 
diffraction breakup is assumed to be similar [S ^| to 
that of stripping. 

The expression (0 gives the contribution of the LM 
distribution coming from each neutron-core state com- 
posing the 15 ground state wave function. For compar- 
ison with the experimental data, we sum up these con- 
tributions weighted by the spectroscopic factors (Table 

The fragment-target interaction cross section is deter- 
mined by the profile function S v {b v ) as © as 



d 2 K(l- \S u {b„)\ 2 ), 



(6) 



where index v denotes the fragment (y = p, 14 N), b v is 
the impact parameter of the i/-th fragment. 

The interaction cross section for the fragmented pro- 
jectile is expressed through the profile functions of the 
fragment-target interaction and the wave function of the 
relative motion of the fragments. 



07 = 



2 J 
1 



MM' J 

J dr** M ,(r) S n S c * JMj (r) 



III. PROFILE FUNCTIONS 

The profile function of the fragment-target interaction 
in Q),©, and JSJ) is determined as an integral of the 
corresponding complex interaction potential 



S v (bv) = exp 



i 

hv 



/ dz V, 



+ z- 



(8) 



where V v t{t) is the fragment-target interaction poten- 
tial, v is the 15 O beam velocity in the laboratory frame. 
The fragment-target interaction potential is determined 
by folding of the fragment density distribution and the 
nucleon-target interaction potential. 

To calculate the nucleon-target interaction potential 
V v x(r) {y — n, p) at energies less than 65 A MeV, we 
use the parameters of the global nucleon-nucleus optical 
potential 01- We also use the interaction potential [f| 
generated from the free nucleon-nucleon (NN) interac- 
tion mm valid at energies from 10 to 2000 A MeV. In 
this case, the nucleon-target interaction potential is ob- 
tained by folding of the target density distribution and 
the nucleon-nucleon interaction potential. For the details 
of the profile function calculations we refer to 0, 01 ■ 

For description of the target and fragment nuclear den- 
sities we use different parametrizations. The 9 Be and 
14 N densities are parameterized in the harmonic oscilla- 
tor model [13 



Pi r ) = Po[l + a(r/a) 2 } exp(-(r/a) 2 ) 



(9) 



The parameter a is related to a |l8| . The parameter a is 
fitted (see the next Section) to reproduce nucleon-nucleus 
and nucleus- nucleus interaction cross sections. 

The 12 C density distribution is approximated by a sum 
of Gaussians 0] as 



p{r) =Y,Ai [e-^~P R ^l^ 2 + e -( r +^«) 2 /7 ; 



(10) 



with the parameters from Ref. |18|. In order to vary 
the calculated cross section obtained with the density 
distribution (|10fl . we introduce a scaling factor (3 and 
replace Ri by (3Ri in (|TU|l . 

All the distributions p are normalized to unity, and po 
is a normalization factor. 

To fit the profile functions, corresponding experimental 
data for interaction (reaction) cross sections on C and Be 
targets at intermediate and high energies are used. 

In the case of 12 C, the scaling parameter (3 in i|10[l is 
fitted to reproduce the experimental data for 12 C+ 12 C 

m m m in m m m m m m an d P + i2 c mm 

interaction cross sections. The best fit is achieved for 
/3=0.94. With this (3 value the 12 C rms radius is 2.37 
fm, that is close to the 12 C rms matter radius 2.33 fm 
obtained in [l9| . 

Figures^, lb show the calculated (dashed gray curves) 
and measured (dots) p+ 12 C and 12 C+ 12 C interaction 
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FIG. 1: The energy dependence of the p+ 12 C and 12 C+ 12 C 
interaction cross sections, 07, calculated with the NN inter- 
acti on p otential. Dots in a) and b) are the experimental data 
[H, [2^ and E^, 113, IH ISa, ISE [2^, [IE E3, Ull , respec- 
tively. The curves correspond to /3=0.94 (dashed gray lines) 
and /3=1 (solid black lines). 



cross sections at energies from 20 to 1000 A MeV. For 
comparison, the cross sections obtained with the charge 
radius of carbon r c — 2.47 fm (/3=1) are also given in 
Fig. (solid black curves). 

The 9 Be density parameter a in © is fitted to re- 
produce the experimental data on the p+ 9 Be [2^, |2|j 
and 9 Be+ 9 Be [23| interaction cross sections. The value 
a = 1.69 fm corresponds the 9 Be rms radius 2.38 fm [l9j . 
To have a measure of sensitivity of the results to the in- 
put parameters of the model, we present the results of 
the calculations with the parameter a — 1.79 fm, also 
allowing a good fit of the 9 Be-nucleus cross section. 

In Figures and|2t> the calculated cross sections are 
compared to the experimental data. These results are 
also compared to the calculations with the Be rms radius 
equal to the Be charge radius, 2.52 fm (a = 1.79 fm) [Isj . 

The p+ 9 Be interaction cross section calculated with 
the NN-interaction potential at energies less than 60 A 
MeV is underestimated, while that obtained with the op- 
tical model potential satisfy the experimental data. At 
higher energies, the cross section calculated with the NN- 
interaction potential is in a good agreement with the ex- 
perimental data |2{jj . 

To test the fitted density parameters of the 12 C and 
9 Be we calculate the interaction cross section in the 



FIG. 2: The energy dependence of the p+ 9 Be and 9 Be+ 9 Be 
interaction cross sections, 07. Dots in a) and b) are the ex- 
perimental data l28l |29| and |23j|. respectively. Dashed gray 
and solid black lines are the calculations with the NN inter- 
action potential with the parameters a — 1.69 and a — 1.79 
fm. The solid grey line are the calculations with the optical 
model potential (OMP). 

9 Bc+ 12 C reaction at the energy 790 A MeV. The value, 
818.7 mb, is very close to the experimental one 806(9) 
mb 0. 

Using the experimental data on the 14 N+p reaction 
[H and the 14 N+ 12 C reaction at the energies 39.3 H3 
and 965 A MeV [l9| we found a = 1.76 fm in the 14 N den- 
sity parametrization 0. This value corresponds to the 
14 N rms matter radius 2.44 fm known from experiment. 

The results of these calculations and the experimental 
data are given in Fig. The cross sections calculated 
with the NN interaction potential are in a better agree- 
ment with the experimental data both for the proton- 
nucleus and the nucleus-nucleus interaction cross sections 
than those obtained with the optical model potential. 

For further calculations of the 15 O break-up on a Be 
target at the energy 56 A MeV we use profile functions 
obtained with the NN interaction potential. 

With the 14 N rms radius we can estimate the 15 O rms 
radius as 

r 2 m C 5 0) = ( r l_ p ) + ^r r 2 „( 14 N), (11) 

where r c _ p is the distance of the valence proton from the 
14 N center of mass, A — A c + A p is the mass number of 
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TABLE II: The calculated (crj) and measured (<jj Xp ) nucleus- 
nucleus interaction cross sections. 



Proj . Target 


E 






° i 




(MeV/ul 

I iVXv> V / Li. 






fmbl 


lb O H Be 


710 




881° 


912(23) 


Be 


710 




920 6 




12 C 


670 




939 


915(13) 


12 c 


710 




945 


922(49) 


a obtained with a = 


1.694 fm 








b obtained with a — 


1.791 fm 








TABLE III: The sin 


gle-particle one-proton removal cross sec- 


tion (crji) and the one-proton removal 


cross 


section (c-p) 


from lo O calculated at the energy 


56 A MeV on a Be target. 


14 N {r,T) 


ot P v 


<T_p 




FWHM 




(mb) 


(mb) 


(MeV/c) 


(l + ,0) 


29.7 


43.3 




178 


(o+i) 


25.9 


10.8 




191 


(l + ,0) 


23.3 


16.2 




198 


(2+,0) 


20.4 


25.6 




209 


Total 




95.9 




191 



the projectile, the valence proton mass number A p = 1. 

With the 14 N rms matter radius r m = 2.44 fm, which 
corresponds (r\ = rjL + 0.8 2 ) to the charge radius 
r c ( 14 N) = 2.57 fm (l8|. and the rms r c _ p distance of 

the proton (r 2 _ p ) 2 = 3.15 fm, the 15 rms matter ra- 
dius is r m ( 15 0) = 2.48 fm. This value is consistent with 
the values obtained in Refs. The corresponding 

15 rms charge radius is r c ( 15 0) = 2.61 fm. 

In Table ITTI the values of the 15 O interaction cross sec- 
tion (0) obtained in the 12 C and 9 Be targets with the 
fitted density parameters are compared to the experi- 
mental data. One can see a good agreement with the 
experimental data [Tflj . 

IV. RESULTS AND DISCUSSION 

The 14 N and 15 O interaction cross sections obtained at 
the energy 56 A MeV on a Be target are ai( 14 N) = 1061 
mb and <7r( 15 0) = 1091 mb, respectively. 

One-proton removal cross sections from 15 O and the 
corresponding FWHM values of the LM distribution of 
the 14 N fragments obtained at the energy 56 A MeV for a 
Be target, are listed in Table ITTT1 All the values are calcu- 
lated with the Be target density parameter a = 1.69 fm. 
The single particle proton removal cross sections, cr_ p , 
and those multiplied by the corresponding spectroscopic 
factors 0]j <7-p, are given for each single particle state. 

The total value of the one- proton removal cross section 
and the LM distribution (last row of Table IlTT|| are found 
as the sum of the proton removal cross sections a^. p and 
the corresponding LM distributions. 

The calculated values of the total proton-removal cross 
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100 1000 



E, MeV/u 



FIG. 3: The energy dependence of the p+ 14 N and 14 N+ 12 C 
interaction cross sections, <tj. Dots in a) and b) are the ex- 
perimental data psl |29| and [TH. |30| . respectively. The cal- 
culations with the NN interaction potential and the optical 
potential (OMP) are shown by solid and dashed lines, respec- 
tively. 




k.., MeV/c 



FIG. 4: Total longitudinal momentum distribution of the N 
fragments (solid line) from the 15 O break-up on Be target at 
the energy 56 A MeV. Dots are the experimental data |4j. 
Dashed line shows the longitudinal momentum distribution 
of the 14 N in the ground state. 
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section and the FWHM (Table |Hl|) obtained with the Be 
target density parameter a — 1.69 fm and the spectro- 
scopic factors are in a very good agreement with the 
experimental values 80 ± 20 mb and 190 ± 10 MeV/c 
Q . With the spectroscopic factors from ^1 the value of 
the proton removal cross section is 92.0 mb. With larger 
target density parameter a — 1.79 fm we get a larger 
values of the cross sections. In this case, the total pro- 
ton removal cross section obtained with the spectroscopic 
factors is cr_ p = 100.2 mb. So one can see that the 
total one proton removal cross section is not very sensi- 
tive to small variations of spectroscopic factors or target 
density parameter. 

In Figure 0] the calculated LM distributions are com- 
pared to the experimental one Note, that the theo- 
retical curves are shifted by 10 MeV/c to the left to be 
compared to the experimental data. 

The solid line in the figure shows the total LM distribu- 
tion obtained with all 14 N states shown in Table ITTT1 The 
dashed line represents the LM distribution from proton 
removal leading to 14 N in the ground state. As it cor- 
responds to the smaller proton separation energy (Table 
ITT}) . the LM distribution is narrower, than that for other 
14 N states. Thus, the FWHM value of the total LM dis- 
tribution is larger than that for the 14 N ground state by 
13 MeV/c. 

The consideration of the 14 N production in the excited 
states does not change significantly the LM distribution 
because each 14 N state (including the ground state) is 
characterized by relatively high proton separation energy, 
and, hence, has nearly the same (Table II 1 1[) LM distri- 
butions. Thus, the value of FWHM is weakly sensitive 
to the weights of the 14 N states and the 14 N excitation. 
However, these contributions are essential in the calcula- 
tions of the proton removal cross section. 

V. CONCLUSION 

In this paper, we present calculations of the one-proton 
removal cross sections from 15 O on a Be target at the en- 



ergy 56 A MeV. The proton removal cross sections, the 
15 O interaction cross section, and the longitudinal mo- 
mentum distribution of the 14 N fragments are obtained 
in the eikonal approximation of the Glauber model with 
the NN interaction potential. In the calculations, the 
production of the 1 N core fragment in the ground and 
excited states is regarded. The calculated FWHM=191 
MeV/c of the total LM distribution is very close to the 
experimentally measured value of 190 ± 10 MeV/c 

The calculated value, 95.9 mb, of the total one-proton 
removal cross section is also very close to the experimen- 
tal value 80 ± 20 mb £| . The break-up cross section is 
about 11% of the 15 interaction cross section. 

Returning to the 17 Ne problem, we see that the con- 
tribution of the proton removal from the 15 O core might 
be essential. In particular, at the energy 66 A MeV (see 
experimental data |2(), we get the cross section of the 
proton removal from the core fragment 94.4 mb. Due 
to the weakly-bound protons blocking the 15 O core in 
17 Ne, this cross section is reduced, contributing about 51 
mb to the total one-proton removal cross section. The 
contribution of the valence proton removal in 17 Ne with 
the spectacular 15 core is about 110 mb Q. Thus, the 
calculated total proton removal cross section will be 161 
mb. This value satisfies the experimental one, 168 ± 17 
mb 2] . Note, that the contribution of the proton removal 
from the 15 O core affects also the width of the total 15 O 
LM distribution. 

As a result, in the reactions with 17 Ne, the proton re- 
moval cross section measured at the energy 66 A MeV on 
a Be target [2j is relatively large compared to the cluster 
model ( 15 0+p + p) predictions Q and the measured 15 
LM distribution is wider than calculated one. 

Therefore, the proton removal from the core should 
necessarily be taken into account in calculations of the 
17 Ne fragmentation. 
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